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Abstract 

We investigate the dynamical behavior of binary fluid systems in two dimen- 
sions using dissipative particle dynamics. We find that following a symmetric 
quench the domain size R(t) grows with time t according to two distinct al- 
gebraic laws R(t) ~ t n : at early times n = 5, while for later times n = |. 
Following an asymmetric quench we observe only n = |, and if momentum 
conservation is violated we see n = w at early times. Bubble simulations con- 
firm the existence of a finite surface tension and the validity of Laplace's law. 
Our results are compared with similar simulations which have been performed 
previously using molecular dynamics, lattice-gas and lattice-Boltzmann au- 
tomata, and Langevin dynamics. We conclude that dissipative particle dy- 
namics is a promising method for simulating fluid properties in such systems. 
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I. INTRODUCTION 



Growth kinetics in binary immiscible fluids has received much attention recently. Phase 
separation in these systems has been simulated using a variety of techniques, including cell 
dynamical systems without hydrodynamics [0] and with Oseen tensor hydrodynamics 0; 
time-dependent Ginzburg-Landau models without hydrodynamics [[|, and with hydrody- 
namics as well as lattice-gas automata P|-|rT|j and the related lattice-Boltzmann tech- 
niques [[T2|,|nj. A central quantity in the study of growth kinetics is the time-dependent 
average domain size R(t). For binary systems in the regime of sharp domain walls, this 
follows algebraic growth laws of the form R(t) ~ t n . In general, for models without hy- 
drodynamic interactions (that is, when there is no conservation of momentum, as is often 
supposed to be the case for binary alloys) the growth exponent is found to be n = |, indepen- 
dent of the spatial dimension. If flow effects are relevant (as is certainly the case for binary 
fluids), and the domain size R is greater than the hydrodynamic length Rh = u 2 /pa, where 
v is the kinematic viscosity, p is the density and a is the surface tension coefficient, then one 
obtains n — | in two spatial dimensions |14]]. In the less commonly observed R < Rh regime, 
two-dimensional lattice-gas automata [11], molecular dynamics ]15[ and Langevin dynamics 
simulations 0,0] find n = |; lattice-Boltzmann studies, by contrast, suggest that n = | flI5| . 
However, the lattice-Boltzmann method does not include any thermal fluctuations which a 
renormalization group approach shows play a crucial role in causing this exponent to assume 
the value of \ ||14|| . In three dimensions, for R < Rh the growth exponent is n = | crossing 
over to n = 1 at late times, with n = | again if R > R h [14| . 

It should be noted that there are still some experimental and theoretical challenges in 
unraveling the behavior of systems in which both the order parameter and the momentum 



are locally conserved ||14|| . Experimentally, for example, it is difficult if not impossible 
to study two-dimensional fluid systems. As far as numerical studies are concerned, it is 
important to recognize that three-dimensional simulations are particularly demanding on 
all the aforementioned techniques and so definitive results are harder to come by than in 
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two dimensions. 

The purpose of the present paper is to take a look at domain growth and phase separation 
in two-dimensional, binary, immiscible fluids using a new simulation technique called dissipa- 
tive particle dynamics (DPD). The basic features of the method are discussed in Sec. 0; here 
we simply note that it is a temporally discrete and spatially continuous, microscopic, par- 
ticulate technique which yields correct hydrodynamical behavior in the macroscopic limit, 
while being easy to extend from two to three dimensions. Comparatively little work has 
been published on DPD and, to our knowledge, nothing at all on its application to phase 
ordering kinetics. We shall find that the method is able to handle domain growth both qual- 
itatively and quantitatively, yielding the correct scaling exponents and displaying a surface 
tension which satisfies Laplace's law. Although the present study is confined to the case of 
two-dimensional systems, we hope to return in the near future with a second paper dealing 
with the three-dimensional case. 



II. WHY DISSIPATIVE PARTICLE DYNAMICS? 

The motivation for the introduction of dissipative particle dynamics by Hoogerbrugge 
and Koelman |16| was to simulate the behavior of complex fluids. Complex fluids include 



fluids in which there are many coexisting length and time scales, and fluids for which a hy- 
drodynamic description is unknown or does not exist at all. Examples are multiphase flows, 
flow in porous media, colloidal suspensions, microemulsions and polymeric fluids. Whereas 
the traditional continuum-based approach to understanding and modeling the hydrodynamic 
behavior of such fluids, based on the formulation and solution of partial differential equa- 
tions, has met with rather limited success, in more recent times new approaches have been 
undertaken, relying on a microscopic description of the fluid in question. In principle, the 
most accurate microscopic approach is based on the use of molecular dynamics (MD) but 
the computational expense of attempting this is so severe that until now only a small num- 
ber of flow phenomena in simple fluids have been achieved, and even then these have been 
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restricted to two dimensions. 

Lattice-gas automata (LGA) have been used as a numerical technique for modeling hy- 
drodynamics since their introduction in 1986 by Frisch, Hasslacher and Pomeau ]T7|] and 
by Wolfram [|T^], who showed that one could simulate the incompressible Navier-Stokes 
equations for a single component fluid using discrete Boolean elements on a triangular lat- 
tice. In essence, LGA dynamics is comprised of two elements: at each discrete timestep, 
particles first collide at vertices on the lattice, the collisional outcomes being controlled by 
local conservation of mass and momentum; then the particles advect freely to neighboring 
sites. Such lattice gas automaton models are computationally much faster than molecular 
dynamics, particularly since the natural time step — the mean free time between collisions — 
is several orders of magnitude greater than that required for MD. The single phase LGA 



method was subsequently generalized by Rothman and Keller |19] to permit the simulation 
of immiscible fluids, and indeed our present work shares certain features in common with 
their model, which has since been investigated with some degree of rigor [p],|10|] . Even more 
recently, the technique has been extended to model amphiphilic fluids comprised of mixtures 
of oil, water and surfactant [ 2D ]. 



However, the LGA method suffers from some disadvantages: the underlying lattice leads 
to the loss of Galilean invariance which, although negligible for creep flows, does present 
problems for flow at finite Reynolds numbers, while the treatment of three-dimensional fluids 
is computationally challenging owing to the complexity of the collisional look-up tables and 
the necessity of deploying a four-dimensional face-centered hypercubic lattice ||. The lack of 
Galilean invariance manifests itself by a spurious factor multiplying the inertial term in the 
momentum-conserving Navier-Stokes equation. For a single-phase lattice gas, this factor can 
easily be scaled away; for compressible flow, or for multiphase flow with interfaces, however, 
the presence of this factor is more difficult to deal with, although various rather involved 



methods have been proposed to remove it pT|,|22 



Dissipative particle dynamics was introduced with the intention of capturing the best 
aspects of MD and LGA; it does away with the problems arising in the latter owing to the 



presence of a lattice, while maintaining the discrete time-stepping element which greatly ac- 
celerates the algorithm compared with MD. Moreover — and this is important from a practical 
perspective — in DPD the extension from two to three dimensions is very straightforward. 
The DPD method involves the motion of massive particles, which are allowed to move in 
space with their positions and momenta described by real numbers. The mass and momen- 
tum of these particles are conserved, but their energy is not. As in LGA, the evolution of 
the model over one time step takes place in two substeps which are continually repeated: (i) 
an infinitesimally short impulse step, and (ii) a propagation step of duration At. Within the 
impulse substep, the momentum of each particle (pj) is modified to reflect its interaction 
with the other particles. During the propagation step each particle coasts with constant 
velocity, completely ignoring every other particle. In mathematical terms, the impulse step 
is described by 



while the propagation step is 



Apj = y^jjejj, (1) 



Aq, = — ( Pi + ApO . (2) 

m 



In these equations m is the mass of each particle, q; denotes the position of particle i, and 
&ij is the unit vector pointing from particle j to particle i. Henceforth, we shall assume 
for simplicity that all particles carry unit mass (m = 1). ^ is a scalar giving the mo- 
mentum transferred from j to i, and in the original model presented by Hoogerbrugge and 
Koelman [IB] has the form 

3 / _ r i3 



n« = { win \ l 'V c ) [n - ^ ■ ^ if r « < rc ' (3) 

if Tij > r c , 

where = |q« — q^l is the distance between particles i and j, and n = N/V is the density 
of the system comprising N particles in a volume V. ( = Uji) is sampled from a uniform 
random distribution with mean and variance Ho. This random component of the momentum 
transfer represents the stochastic effect of the collisions and gives rise to a fluid pressure, 



while the second, dissipative, term inside the square brackets of Eqn. (|3]) is responsible for 
the fluid viscosity. The temperature is controlled by H^, whose variance is a measure of 
the thermal fluctuations in the system. Note that r c , which occurs in the factor multiplying 
the one in square brackets in Eqn. ([$]), is a cut-off radius beyond which no interaction 
(momentum transfer) is possible. The presence of this cut-off makes the interactions local 
and the DPD algorithm correspondingly fast. 



The property of detailed balance is satisfied by DPD [^3| and so a Gibbsian equilibrium 



state is guaranteed to exist. However, the statistical mechanical analysis of Espanol and 



Warren p4] confirmed that the original DPD model presented by Hoogerbrugge and Koel- 
man |16| does not lead to the physically correct equilibrium distribution, a feature which is 
particularly marked when the relative time step At > 1. 

This choice of time step (At = 1) has, nevertheless, been used in almost all of the 
DPD work so far reported on |16lE5]-|27j . Two simple modifications to the basic model are 



suggested by Espanol and Warren |H[ which ensure that the DPD equilibrium state is the 
canonical ensemble. The first modification is reducing the time step length (a factor of 
ten is sufficient), and the second is the inclusion of an extra factor of 2(1 — ry/r c ) into the 
dissipative term in the change of momentum Eqn. (|3|), so that the momentum transfer scalar 



ttij becomes 



n 



irrin 



Uij - 2u ( 1 - ) (pi - p i; ^, 



(4) 



This alteration guarantees that the model obeys a fluctuation-dissipation theorem which is 
very similar to that obtained in conventional Brownian motion |p4|| . The theorem enables 
one to relate the amplitude of the noise (that is, the fluctuations in 11^) to the temperature 
of the system. 

In the original paper by Hoogerbrugge and Koelman, it was stated (but not demon- 
strated) that the DPD model of a simple one-component fluid satisfies the Navier-Stokes 
equations in the mean-field limit [HJ. Espanol has recently explicitly derived the hydrody- 



namic equations for the mass and momentum density fields in DPD ||23||. However, these 



equations are not the central results of Espanol's paper, since the DPD properties of mass 
and momentum conservation, coupled with Galilean invariance and the isotropy of the mi- 
croscopic equations of motion, effectively guarantee that at a macroscopic level the Navier- 
Stokes equations will emerge. What is more significant is the correction of the original 
simple expressions for the speed of sound and the kinematic viscosity to more complicated 
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III. A DPD MODEL FOR BINARY IMMISCIBLE FLUIDS 

Immiscible fluid mixtures exist because individual molecules attract similar and repel 
dissimilar molecules. The most common example of such behavior arises in mixtures of oil 
and water; the non-polar, hydrophobic molecules of oil attract one another through short 
range van der Waals forces, while the polar water molecules enjoy more complex, long-range 
hydrophilic attractions which are dominated by electrostatic interactions including hydrogen 
bonds. At the atomistic level employed in molecular dynamics, such interactions demand a 
detailed treatment. However, to obtain accurate meso- and macroscopic level descriptions 
using DPD, the microscopic model can be drastically simplified. To model the interactions 
of dissimilar particles in a binary immiscible fluid within DPD, the simplest modification 
to the single phase DPD algorithm is to introduce a new variable, called the "color" (by 
analogy with the Rothman-Keller model), which has two possible values — for example, "red" 
for oil and "blue" for water. Identical interactions are used between particles of the same 
color, while we increase the mean and variance of the random variable IT^ when particles of 
different color interact, which has the effect of creating a repulsion between particles from 
the two different phases p5| . That is, 



U[0, 2IIo] if particles are of the same phase, 

% e { (5) 

U[0, 2(II + n rep )] if particles are of different phases. 
It should be emphasized that Eqn. (|5]) is a minimal modification to the single-phase 
DPD model and is symmetric under the interchange of particle colors. It would be quite 
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straightforward to generalize this to the asymmetric case by, for example, making the mean 
n of the stochastic terms different for the two colors and/or likewise adjusting the dissipative 
terms through the selection of different values of to in Eqn. (^). In the present paper, 
however, we shall not consider such situations, preferring to concentrate on the properties 
of the simpler model implied by Eqn. (^j). It has been shown that detailed balance is also 
satisfied by the two-phase DPD model |[28|| . As for the single- phase DPD fluid, one can show 



that, macroscopically, the Navier-Stokes equations are obeyed within homogeneous regions 
of each of the two immiscible fluids. 

The implementation of the DPD algorithm is very similar to that of conventional MD 
algorithms ||29|| . For example, the periodic spatial domain (the simulation cell) is divided 



into a regular array of equally-sized link cells, such that each side of the rectangular domain 
has an integer number of cells and each cell is at least r c across. Each link cell consists of 
a dynamically-allocated array of particles and pointers to the neighboring cells. Individual 
particles consist of the position-momentum vector pair and a color index. 

For each time step we calculate first the impulse and then the propagation step, as 
described in Eqns. (|IJ) and (0). In the impulse step we iterate through the particles in 
each link cell, calculating the change in momentum of each particle as it interacts with 
the particles in the same and neighboring link cells. Since the momentum is modified by 
particle pairs we need to ignore half of the neighboring cells to avoid duplication. When 
considering a new particle pair we first compare the square of the separation distance with 
rf, and skip to the next particles if the pair is out of range. In the propagation step we 
iterate through the particles in each link cell again, allowing them to coast for time At. The 
complete state of the system may be written to file, and other calculations to determine for 
example the temperature and pressure of the system can also be performed. Given constant 
r c and number density n = N/V the system scales linearly (in both computation time and 
memory size) with increasing number of particles, N. To give an example, the calculation 
of the motion of a system of 40,000 DPD particles with number density n = 4 for 10,000 
time steps takes 11 hours CPU time and 13 MB memory on an 133 Mflop/s (theoretical) 
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DEC Alpha. 



IV. RESULTS 
A. Scaling Laws for Binary Fluid Separation 

Several simulations were run to model the binary separation of a mixture of two immis- 
cible phases. By definition, for studies of symmetric quenches, exactly half of the particles 
were in the "red" phase and half in the "blue" phase. Asymmetric quenches were studied 
with the ratio of the numbers of colored particles at both 60:40 and 70:30. The initial state 
of the system was completely random; the positions of red and blue particles were chosen 
from a uniform distribution, and their velocities were chosen from a uniform distribution in 
direction and a Gaussian distribution in magnitude. Physically, this initial state corresponds 
to starting with the system quenched in temperature from a state above the spinodal point 
at which the fluids are miscible. 

A reasonably large number of particles (40,000) was used in order to enhance the sta- 
tistical accuracy of the data obtained during these simulations. Care was taken to ensure 
the box dimensions were suitable for accurately simulating the phase segregation process 
for large times without the size of the domains becoming close to that of the periodic box, 
thereby introducing artifacts into the observed behavior. 

Simulations started from a symmetric quench were allowed to evolve for 10,000 time steps 
and the asymmetric quenches were evolved for 100,000 time steps. The state of the system 
was recorded at regular intervals throughout. Fig. [I] shows the state of a single system at 
six different snapshot times following a symmetric quench. At each time for which the state 
of the system was recorded, the structure function 



S(k,t) 



i J (p 6 (x, t) - p r (x, t) - (p b ) + (p r )) exp (— 2ttz k • x)dx 



(6) 



was calculated. In this equation the subscripts b and r indicate blue and red particles, 
so that p&(x, t) is the (spatial) mass density of the blue particles at time t, and (pb) is 
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the average mass density of blue particles. The use of the structure function to characterize 



binary fluid phase separation is widespread [H11.3C . Note that, because of the imposition of 



two-dimensional periodic boundary conditions on the simulation cell, the structure function 
is only meaningful when evaluated at points in the Fourier space satisfying 



f m n \ 
k = (k x ,k y ) = , — J ; m,neZ. 



(7) 



where L x and L y are the lengths of the box in the x- and ^/-directions of real space, respec- 
tively. 

To extract the quantity best characterizing the time-dependent domain size of the state 
of the system, the mean of k weighted by the structure function is often calculated. However, 
unlike the case of lattice-based simulation techniques, the positions of the particles in DPD 
are described by continuous variables rather than discrete points on a grid, and so the 
structure function is meaningful even as |k| approaches infinity. The nature of the system 
is such that the structure function has a finite value throughout k-space, approaching a 
constant value at large distances from the origin, so that the mean is not defined. Since 
we are not interested in the asymptotic value of the structure function, a function -F(k, t), 
which vanishes far from the origin, was fit to the grid of points containing the meaningful 
values of S(k,t). The function that we chose as the best fit for S(k,t) in these simulations 
is 

F(k,t) = F(\k\,t) = c Q (t) |k|exp (- Cl {t) |k| 2 ) . (8) 

Note that F(k, t) is taken to have radial symmetry as we are not interested in the orientation 
of the domains in real space. A sample structure function with its fit F(h,t) are shown in 
Fig. |[ The reciprocal of the mean of |k| weighted by F(k, t) is denoted R(t), and is 
interpreted as the domain size characterizing the state of the system at time t. 

By plotting our computed values of R(t) on a log-log plot versus t, it is easy to see any 
change in the exponent of the scaling law 

R{t) oc t n . (9) 
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To obtain accurate results it is important to ensemble-average several simulations differing 
in their random initial conditions. For each simulation, the log-log plot was examined 
for possible exponents. From each of the symmetric quenches it was possible to mark 
the existence of two separate scaling regimes, with the crossover at roughly t = lOOOAt. 
Only a single scaling regime was observed from the asymmetric quenches, even though 
the simulations were allowed to evolve an order of magnitude longer than the symmetric 
quenches and the final domains were significantly larger. For each regime in all simulations 
a straight line was fit to the data to determine the exponent. The complete results for 
symmetric quenches are listed in Table |, while a sample log-log plot is shown in Fig. [3[ All 
of the figures in this section are derived from one and the same simulation, and so may be 
directly compared. The results in Table | strongly suggest a crossover from n = \ to n = |, 



which is in agreement with results from lattice gas automata | IT| , molecular dynamics fl5 



and Langevin dynamics simulations |7j as well as a renormalization group analysis JT4|. A 
growth exponent of n = ^ was observed from the asymmetric quenches, independent of 
the ratio of particle numbers. This exponent was observed at all times in these off-critical 
quenches, just as in previously reported Langevin simulations H0]. 

Some short simulations (3,000 time steps long) were also run to determine the early- 
time exponent for a symmetric quench where momentum is not conserved. The momentum 
transferred between a pair of interacting particles was violated by adding a vector 15-20% 
of the magnitude of the interaction and random in direction (chosen from a uniform distri- 
bution) to each particle's momentum. The growth exponent observed in these simulations 
was n = 0.329 ± 0.005, in agreement with theory |L4 . 



The complete set of model parameters used in these simulations is listed in Table 0. In 
this table, L x and L y are the box dimensions, while n is the number density N/V, where 
V is the two-dimensional volume of the box. The value of n = 5.0 was chosen since the 
immiscible phases then regularly reach complete phase separation asymptotically and yet 
the system does not behave too randomly. 
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B. Bubble Surface Tension 



The detailed way in which phase separation occurs in binary fluids depends, among 
other things, on the interfacial tension which exists between the two immiscible phases. In 
particular, as noted earlier, when hydrodynamic effects are important, the crossover between 
diffusive (Lifshitz-Slyozov) and hydrodynamic regimes should occur when the domain size 
R(t) is larger than the hydrodynamic length Rh = v 2 j per, where v is the kinematic viscosity, 
p is the fluid density and a is the surface tension coefficient. A further important test of our 
DPD model for binary fluid separation is thus to check on the existence of a surface tension 
between the two phases by confirming the validity of Laplace's law using a series of bubble 
simulations inspired by earlier lattice gas analogues [fTT|j3lJl . 



As with the domain growth simulations of Sec. [IV A| , 40,000 particles were placed in 
a two-dimensional periodic box of the same dimensions as those listed in Table |J, but 
the simulations were now run for 40,000 time steps. (For a few of the smaller bubbles, 
simulations were run with only 6,400 particles and to 5,000 time steps. For these smaller 
bubbles accurate results could be obtained without the additional computation required 
by the larger system.) The results here have the same model parameters as in Sec. [IV A| 
(see Table [TT]) with the sole exception of n , which was set to 0.25 in order to increase the 
signal-to-noise ratio. The initial state of the system has all the particles placed with random 
positions and velocities, but the particles are red within a radius Rq of the center of the box, 
and blue outside this distance. This bubble changes rapidly within the first few time steps, 
but settles down to a state approximating equilibrium after about 8,000 steps for the larger 
simulations (2,000 time steps for the smaller simulations). Fig. ^ shows one such bubble at 
equilibrium. 

These bubble experiments thus enable us to compute the equilibrium value of the inter- 
facial tension between red and blue phases. Equilibrium statistical mechanics tells us that 
we can write the instantaneous pressure function (V) of a system in terms of the internal 
virial (W) and the instantaneous temperature function (T) p9| as 
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V = ^^, (10) 

where 

Nk B T=\j2—> ( n ) 

3 mj 

and 

where /cb is the Boltzmann constant. The thermodynamic pressure (P) and temperature 

(T) of the system are the time averages of V and T, respectively It is possible to compute 



a temperature because the important property of detailed balance is satisfied by DPD [23 
the fluctuation dissipation theorem relates this to the noise. (We note in passing that 
Hoogerbrugge and Koelman [Tj| gave an equation of state for a homogeneous DPD fluid 
which relates the fluid pressure to the fluid density. However, since the theoretical basis for 
this equation is not adequately explained we have preferred to calculate the pressure using 
statistical mechanical first principles.) 

By considering only the blue particles far from the interface (r > 1.3i?o) we can calculate 
the thermodynamic pressure in homogeneous regions outside the bubble (P ou t)- Similarly, 
we can calculate the pressure inside the bubble (Pj n ) by using particles with r < 0.7i?o- For 
a circular bubble in two dimensions, Laplace's law says 

Pin Pout = ~~f^~- (13) 

To verify this we calculate the mean pressure difference at several values of Ro, and plot 
the pressure difference Pi n — P out versus the reciprocal of Ro. For each of the simulations 
the equilibrium pressure difference was averaged to a single mean value. The mean pressure 
difference and corresponding error for each bubble size reported in Fig. |5| are the mean and 
standard error of the mean for the set of simulations at each bubble size. We can see from 
our results in Fig. |5] that we have the expected linear behavior. The slope of this line is the 
surface tension coefficient, a. 
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V. CONCLUSIONS 



We have studied binary immiscible fluid behavior in two spatial dimensions using a 
novel simulation technique called dissipative particle dynamics. We have found algebraic 
scaling laws in agreement with expectations [T%\, the two-dimensional growth exponents 
being | and | at early and late times respectively for symmetric quenches, and \ throughout 
for asymmetric quenches. Symmetric quenches for which momentum transfer was violated 
displayed an early-time growth exponent of |. This scaling behavior has previously been 
observed in molecular dynamics [I3|, Langevin dynamics and lattice gas automata |TT 



simulations, and is also in accord with the results of a renormalizat ion-group approach which 
takes into account the mechanism of droplet coalescence due to noise-induced Brownian 



motion | l3fl . We have also verified that Laplace's law holds in a series of simple bubble 
experiments, confirming the existence of a surface tension between the two phases. 

We conclude that dissipative particle dynamics is a promising method for simulating the 
properties of fluid systems. 
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FIGURES 



FIG. 1. Binary Separation Particle Positions 
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FIG. 2. Structure Function and Sample Fit 
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FIG. 3. Sample R vs. t log-log Plot 




FIG. 5. Pressure Difference vs. 1/R 
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TABLES 



TABLE I. Binary Separation Exponents for Symmetric Quenches 


Simulation 


Slope 1 


Slope 2 


logt 


logR 


number 






crossover 


crossover 


8 


0.4701 


0.648 


3.000 


0.865 


9 


0.4648 


0.640 


3.000 


0.869 


11 


0.4805 


0.639 


3.000 


0.886 


12 


0.4788 


0.642 


3.114 


0.932 


13 


0.4804 


0.649 


3.079 


0.916 


14 


0.4859 


0.627 


3.000 


0.884 


15 


0.4759 


0.676 


3.114 


0.929 


16 


0.4717 


0.673 


3.079 


0.912 


17 


0.4855 


0.696 


3.114 


0.943 


Mean 


0.477 


0.65 


3.06 


0.90 


StdDev 


0.007 


0.02 


0.05 


0.03 
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TABLE II. Binary Separation Model Parameters 



Model 




Parameter 


Value 


At 


0.1 


m 


1.0 


n 


4.0 


AT 
iV 


4UUUU 




2.0 


n 


5.0 


TT 

li- re p 


U.6 


r c 


1.3 




100.0 


Ly 


100.0 
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